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Abstract 

In this paper we evaluate finite temperature corrections to the dark matter 
relic density within the context of minimal supersymmetry with a neutralino 
LSP. We identify several regions of parameter space where the WIMP annihila- 
tion cross section is especially sensitive to small corrections to the undelrlying 
parameters. In these regions, finite temperature effects have the potential to 
be important. However, we shall show by explicit calculation that these ef- 
fects are small. In the regions we investigated, the maximal corrections are on 
the order of 1CP 4 and are therefore negligible compared with theoretical and 
experimental uncertainties. 
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1 Introduction 



It is now well accepted that approximately 20% of the matter in the universe is 
in the form of a new non-baryonic neutral and weakly interacting particle whose 
microscopic properties have yet to be determined. This material has been termed 
dark matter. Experimentally, the ratio of the dark matter density, on cosmological 
scales, to the critical denstiy has been measured within 10% to be Qdm ~ 0.2 [T|. 
Given a theory of dark matter, the relic density can be predicted by caculating the 
particle's annihilation cross sections and integrating the Boltzmann equation through 
decoupling to the present day. It is interesting to note that if we assume that the dark 
matter particle annihilates through weak interactions and has a mass on the order 
of the electroweak scale, we obtain the correct value for the relic density (see, for 
example, |2j). This fact has led many theorists to concentrate on weakly interacting 
massive particle (WIMP) models of dark matter. 

Assuming that dark matter is composed of WIMPs that are in thermal equilibrium 
in the early universe, their density today is given approximately by the expression |3] 



where Sq is the current entropy density of the universe, p c is the critical density, 
~ 80 is the number of relativistic degrees of freedom at the time of freeze-out, (av) 
is the thermally averaged annihilation cross section and Xf = m/T is the ratio of 
the WIMP mass to the freeze-out temperature. For weak scale interactions Xf ~ 25. 
From the above value of Xf it is clear that the WIMP is non relativistic at the time 
of freeze-out and for this reason it has long been supposed that finite temperature 
effects can be neglected when calculating the annihilation cross section. Indeed, these 
corrections tend to be suppressed either exponentially, or by several powers of T/m 
(see, for example, j3j). However, currently favored models of dark matter can often 
lie in regions of parameter space where the annihilation cross section is extremely 
sensitive to small variations in the underlying parameters. For this reason we consider 
the issue of finite temperature corrections to the relic density worth revisting. 

For definiteness we consider models of dark matter within the context of the min- 
imal supersymmetric standard model (MSSM). Dark matter scenarios in the MSSM 
have been discussed by many authors. In paricular, the parameter space of a sub- 
class of models - minimal supergravity (mSUGRA) - has been extensively mapped 
out [7J |H] • Viable mSUGRA models fall roughly into four parameter regions: the 
bulk region, the focus point region, the coannihilation region and the A-funnel. In the 
bulk region the superpartners are all relatively light and the annihilation is dominated 
by lepton production through P-wave t-channel and s-channel slepton exchange. In 
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the focus point region the scalar superpartners are very heavy. In this case the LSP 
obtains a significant higsino and wino component and the annihilation is dominated 
by s-channel W and Z production. The coannihilation region is characterized by a 
nearly degenerate next-to-lightest super partner (NLSP). In this scenario the NLSP 
has a finite density during freeze-out and coannihilates with the LSP. Finally, the 
A-funnel region is characterized by ~ 2771^ where tua is the mass of the CP-odd 
heavy Higgs A and is the mass of the LSP. In this region of parameter space 
the WIMP annihilations proceed mainly through resonant s-channel A exchange. 

Three of the above regions - the bulk, the coannihilation region and the A funnel 
- exhibit sensitivity to one or more of the underlying MSSM parameters. In the bulk, 
the S-wave annihilation is helicity suppressed. The dominant channel is thus P-wave, 
which is itself suppressed by two powers of the WIMP velocity, v. Any enhancememt 
to the helicity breaking masses of the standard model fermions would thus significantly 
enhance the overall annihilation cross section. In the A-funnel the WIMP annihilation 
cross section is naturally highly sensitive to the width of the A boson. Finally, in 
the coannihilation region the cross section is sensitive to the abundance of the NLSP, 
which is exponentially dependent on it mass. In all these cases small corrections due 
to finite temperature effects have the potential to be significant. 

Similar issues, in which dark matter annihilation cross sections are tuned by the 
adjustment of mass differences or particle mixings occur in many other models, both 
within and outside of the context of supersymmetry. 

In this paper we will present specific estimates of the finite temperature effects on 
the relic denstiy in the scenarios that we have discussed. The outline of the paper is 
as follows: In Section |21 we review elements of finite temperature field theory relevant 
to our analysis. In Section |3] we calculate the thermal corrections to the abundance 
of the NLSP at freeze-out in the stau coannihilation region. In Section 0] we analyze 
the enhancements to the helicity breaking standard model fermion masses and their 
effect on the relic density. In Section we evaluate the thermal corrections to the 
width of the A boson in the A-funnel region. Finally, in Section El we review our 
analysis and present conclusions. 

Other corrections to relic density predictions stemming from a better understand- 
ing of the high temperture environment in the early universe have been suggested. 
Most recently, Hindmarsh and Philipsen jB] have analyzed the deviation of the QCD 
equation of state from that of an ideal gas. They found that this effect could lead to 
measurable a correction to the dark matter abundance. It is, however, beyond the 
scope of this paper to discuss these recent developments. 
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2 Thermal Field Theory 



At finite temperatures, field theory calculations must reflect the fact that all correla- 
tion functions are thermally averaged. The expression 



for the zero temperature vaccum expectation value of an operator must be replaced 



where the sum is over a complete set of states. Perturbation theory then leads to a 
set of corrected Feynman rules for S-matrix elements. To a good level of accuracy 
the corrections to relic density calculations can be understood through the thermal 
shift in physical parameters such as the masses and widths of interacting particles. 
For this reason we will concentrate on the corrections to the real and imaginary parts 
of propagators and we do not find it necessary to evaluate general finite temperature 
Green's functions. 

Two parallel formalisms exist for perturbative calculations at finite temperatures; 
the imaginary time formalism (ITF) which was first introduced by Matsubara jU] 
incorporates a contour along the imaginary time axis to obtain the thermal average. 
This technique is best suited for the calculation of equilibrium properties of hot 
plasmas but can be extended through analytic continuation to real time correlation 
functions. The real time formalism (RTF) was developed later by Schwinger |lUj . 
Mills [TT] and Keldysh jT^j and incorporates dynamics more naturally. The latter 
also provides more calculational ease. It is this technique which we shall utilize in 
the following. 

In the RTF the degrees of freedom in a theory are doubled. The new fields do not 
create or annihilate external particles but may appear in loops. The new fields do 
not couple to the physical fields through Feynman diagram vertices but rather mix 
with them through propagators. It follows that these new fields do not appear at one 
loop level, and we can ignore them in this analysis. Then, the propagators of scalars, 
chiral fermions, and gauge bosons, respectively are [T3] 
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We have used Feynman gauge and the fermionic propagator is given in two com- 
ponent notation. The above equations contain the Bose-Einstein and Fermi-Dirac 
distributions, respectively, 

1 



n(p ) 
n(p ) 



e P\po\ - 1 
1 



Throughout this paper, upper case letters will denote four momenta and lower case 
letters will be reserved for the spatial components. If the fermionic Lagrangian con- 
tains a helicity breaking mass term we will obtain a corresponding Feynman rule for 
a mass insertion, given by 



- 2tt5(P 2 - m 2 )n(p 



(5) 



P 2 + irj 

where c is the antisymmetric 2x2 matrix defined so that c±2 — — 1. To evaluate the 
imaginary part of self energy diagrams we will use a generalization of the zero tem- 
perature cutting rules, due to Kobes and Semenoff [Hj. As in the zero temperature 
case, the imaginary part of a diagram is obtained by replacing full propagators by 
"cut" propagators. For fermions these are given by 

S>(P) = 27re(p )(l-f(p ))6(P 2 -m 2 ) (6) 
S<(P) = -2ne(p )f(p )5(P 2 -m 2 ) (7) 

where e is the sign function and 

f ~ {po) = ^T+T 

Note the absence of the absolute value in the argument of the exponent. 

Throughout this paper, we will take the plasma to be comprised solely of light 
(rrif < 10 GeV) standard model fermions and massless gauge bosons. We will consis- 
tently neglect terms representing interactions with superparticles, heavy gauge bosons 
and heavy Higgs bosons in the thermal plasma. This is a reasonable approximation 
since such terms will be exponentially suppressed by the ratio of their mass to the 
freeze-out temperature. 



3 The LSP-NLSP Mass Splitting 

An often considered scenario for coannihilation involves the stau as the NLSP. In this 
section we will analyze the finite temperature corrections to the model. The anni- 
hilation rate is exponentially dependent on the mass splitting between the stau and 
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Figure 1: One loop diagrams contributing to the self energy of the stau. 

the lightest superpartner. Thus it is interesting to compute the thermal corrections 
to this mass splitting. 

The full sfermion propagator is given by 



K 2 -ml- U(Ky 



Where H(K) is the one paricle irreducible self energy The particle's dispersion rela- 
tion, ko(k) is obtained by solving the equation 

kl - k 2 - m 2 - Yl(K f = 

and the mass shift is then given by the limit 6m = k^ — mo . At one loop level, the 
thermal corrections to the stau mass originate from the Feynman diagrams depicted 
in Fig. ^ Each propagator has two terms, shown in eqs. ©-Q. The first represents 
a contribution from a virtual particle and the second from a thermal particle. The 
first term is the zero-temperature component of the propagator. To obtain the ther- 
mal contributions to the real part of a one loop diagram we need to consider terms 
involving one thermal and one virtual field. As mentioned, terms involving thermal 
heavy particles will be suppresses by small statistical factors and we can safely neglect 
them. 

We begin with the photon loop. Using (J2J) and (j3J), 

_ ilp -ioop = (ie)2 J AL_(K -P + KYD^{P)D{P -K){K-P + K) v (8) 
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We write the real part of the self energy and ignore the T = terms as well as those 
representing interactions with sfermions in the plasma. This gives 



ReTPf loop 



ie 2 Q 2 



d 4 P 

(27T) 4 



(AK 2 + P 2 -4K -P) 
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where m is the zero temperature mass of the stau and the subscript (3 denotes the 
thermal part of the diagram. This can be evaluated to give 



ReIff loop 



~vJ dvv 



2 + (Lf(m ) + L^(m )) 



K 2 + ml 
Akp 



n{p) (10) 



where, following ^H], we have defined 



Lf(K,p,m) = ±ln 



Note that 



lim — — 



'K 2 -m 2 ±2(k + k)p\ 
K 2 — m 2 ± 2(k - k)p) ' 

4p 



11) 
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k-^o k ko ± 2k p — ml ' 

To first order in a we can set ko(0) = m inside the integral, which then considerably 
simplifies. We obtain for the mass shift 



8{m 



2\7~ loop 



-- / dppn(p) = ~^-T 2 . 
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The expression for the photon tadpole is 



j^j-jtadpole 



d 4 p 
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Evaluating this and taking the limit as k — > we obtain 
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Next, we compute the lepton-gaugino loops. Using (jSJ) and following the same 
steps as those for the photon loops we obtain for a digaram containing virtual neu- 
tralinos 



5(m 



2\tau-N 
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(16) 



(1-fi 2 ) 2 - Ay 2 e»/* + l 
where G t n 1 are the couplings of the tau stau and neutralino, ^ = mi/m? is the mass 



of the i neutralino scaled by the stau mass, /i T = m T /m f , y = Jx 2 + fi 2 , t = T/mf 
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Figure 2: Feynman diagram contributing to the neutralino self energy at one loop 
order and in the limit of a plasma populated by light SM particles. 

and the sum runs over the four neutralino types. In the derivation of this expression 
we have consistently neglected terms that are proportional to the neutralino statistical 
factor. A similar expression gives the contribution from the diagrams with virtual 
charginos. 

In the limit of the plasma being populated solely by light SM particles, the neu- 
tralino obtains a thermal mass shift from self energy diagrams of the type depicted 
in Fig. with a virtual sfermion and a thermal standard model fermion. Note that 
we have neglected helicity flipping diagrams since these will be suppressed by light 
fermion masses. 

-*£ = -(*e) 2 X> /JV | 2 / 0^ 2 cS(P)cD(K + P) (17) 

where the sum runs over all fermion species inhabiting the plasma and Gf^ is the 
neutralino- fermion-sfermion coupling. 

In order to extract the mass shift we should better understand the structure of 
this expression. The most general form of the 4-component self energy allowed by 
the Majorana nature of the neutralino is 

U(K) = aKrf + bu^l" + c (18) 

where is the 4-velocity of the thermal bath. In the approximation that we ignore 
the helicity flip diagrams, c = 0. Then, the dispersion relation, 

det(K ■ 7 - E 4 (iT) - m ) = 

taken at zero momentum, leads to a thermal mass given by the equation 

(1 -afm 2 = {m Q ) 2 + b 2 . (19) 

Examining (fT7|) we see the Tt(K~) = — £(— K) so we can conclude that b = 0. The 
identity 

5m = am = —Tr {Re E} . 
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then follows. Evaluating this, we obtain 



TT 



2am N 2 r x z y 1 
Sm = — \ G Jn\ dx- -j-, — - . (20) 



(l-/i 2 ~) 2 - V e»/* + l' 



where here y = ^Jx 2 + fi j. 

Using ([EH) , P^fy an d ff2T?j) We can calculate the correction to the stau-neutralino 
mass splitting. A typical, well studied, point in the coannihilation region is given 
by benchmark point LCC3 which is described in detail in [T7| . At this point the 
mass splitting between the stau and the lightest neutralino is approximately 10 GeV. 
In Fig. El we plot the total shift in this mass splitting, at LCC3, as a function of 
temperature. At freeze-out, which occurs at approximately 5 GeV, thermal effects 
account for a fractional correction of about 10~ 4 . The WIMP annihilation rate is 
determined by the relative Boltzman factor between the two particles, which is then 
corrected by 



e 



-(m f -m N )/T _^ e -(m f -m N )/T e -(m f -m N )V) 4 /T ^ e -(m f -m N )/T ^ — 2 X 10~ 4 ) (21' 



We evaluated the effect of the shift in mass splitting on the relic density using the 
DarkSUSY [THj package. We found, as expected, that the change in the relic density 
prediction is roughly 

, n 5 (Am) 

— f» 1.0 x -. (22) 

U DM Am 



4 Helicity Breaking Fermion Masses 

It would seem that a straightforward way for two WIMPs to annihilate would be 
by production of fermion pairs through t-channel and u-channel sfermion exchange. 
For gravity mediated supersymmetry breaking this is especially true since the LSP 
is typically almost pure bino and as such does not couple to the Higgs boson or to 
gauge bosons. But S-wave annihilation of neutralinos into fermion- ant if ermion pairs 
is helicity suppressed by the square of the ratio of the mass of the fermion to the mass 
of the LSP PHI- In the bulk region of the mSUGRA parameter space, this number 
is on the order of 10~ 4 for tau or bottom quark production. P-wave annihilation is 
not helicity suppressed but is proportional to two powers of the WIMP velocity. This 
is typically about (T/ttiat) w 10~ 2 . Thus, the helicity suppressed S-wave processes 
contribute to the annihilation at a percent level. 

Thermal effects cannot on their own give fermions helicity breaking masses but 
they do have the potential to enhance (or suppress) existing mass operators. The 
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Figure 3: Thermal corrections to the masses of the stau (solid) and lightest neu- 
tralino (dashed) in the stau coannihilation region. The values are scaled by the 
zero-temperature mass splitting. 



thermal corrections to these operators must vanish as the zero temperature mass 
becomes small and will thus be proportional to it. It is possible that the coefficient 
of proportionality could be large since the SM fermions are light compared to the 
ambient temperature. 

Assuming that only the light gauge bosons and SM fermions are thermal, the 
diagrams contributing to the fermion helicity breaking masses are those depicted in 
Fig. EJ For the photon loop 



<S = (*e) 2 / S^i^S(K - P)D{P), v a\ 



(27T) 



(23) 



We evaluate the integral and take the real part of E, neglecting the zero temperature 
terms, to obtain 



ReZ t 



a 



7T 



-m/c / d 4 P 



2it5(P 2 



2ti5(P 2 - m 2 



n(p ) (24) 



(K + P) 2 -m 2 (K + P) 2 
where is the fermion mass. We can perform the po integration using the delta 
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Figure 4: Feynman diagrams contributing, at one loop order, to the helicity breaking 
mass of the standard model fermions. 



function. The first term in brackets yields 
Re^P = ^m f cf d3p 



P 7T J p 



+ TP?. na„ i rzzr^z —2 n (p) 



K 2 + 2(k — k cos 6)p — mj K 2 — 2(k + k cos 0)p — mj 

Oi f 

= —rrrifC I dp(Lf(m f ) + L{ (m /))n(p) 
irk J 

(25) 

Since the fermions produced in the WIMP annihilation are on-shell, we may use 
K 2 = m 2 inside the integrals. It is then clear from (fTT]) that the integrand vanishes. 
So, we see that at one loop order massless gauge bosons cannot enhance the fermion 
mass operator. As we will see, all finite coorections will be due to massive particles 
populating the plasma. This brings us to the second term in (|24j) which yields a similar 
expression that is non zero only for fermions with a non-negligible mass. Thus, we 
consider only tau and bottom quark production. For the tau 

Re^f = i^rn f c J d P ^(L+(m f ) + L^(m f ))n(E). (26) 

where 

+ / s (K 2 + m 2 f ±2(k E + kp)\ , s 

Lt(K,p,E,mf) = ±\n[— { ^ . 27 

and E = ^Jp 2 + mj. This integral should be evaluated for on-shell fermions with 
momentum k on the order of the LSP mass. In the above expression, the coefficient 
of mfC represents the correction to the mass operator. This term, too, would vanish 
identically but for the effect of the non- zero fermion masses. For the tau, this is the full 
result. The bottom quark receives another, stronger correction from thermal gluon 
loops. These will be equal to the above expression with an addional multiplicative 
factor of (4/9)(a s /a). For the purpose of this calculations we have considered the 
LCC1 benchmark point [201 as a cannonical point in the bulk. We have evaluated 
8m/ m numerically for the bottom and tau. The results are plotted in Fig. El We see 
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Figure 5: Thermal correction to the helicity breaking mass operator of the bottom 
quark (solid) and tau lepton (dashed). 

that, at freeze-out, the corrections to the fermion mass operators are on the order 
of 10~ 4 and are thus negligible compared to existing theoretical and experimental 
uncertainties. 



5 The Width of the A Boson 

In the so called A-funnel region of mSUGRA parameter space, the mass of the CP- 
odd heavy Higgs boson (A ) is roughly equal to twice the mass of the LSP and the 
WIMP annihilation cross section is dominated by resonant s-channel annihilation to 
bb, t + t~ . The value of the cross section is, therefore, highly sensitive to the width 
of the A . We would like to investigate the thermal corrections to this quantity and 
their effect on the relic density predictions. 

In the limit that the plasma is populated only by massless gauge bosons and SM 
fermions, the A self energy receives corrections from a single Feynman diagram - the 
fermion loop depicted in Fig. |Hl We would like to evaluate the imaginary part of this 
diagram. As discussed in Sec. El we use the finite temperature cutting rules from [T3] 
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Figure 6: Feynman diagram contributing to the A boson self energy. 

to do so. Following the notation in [2T], we calculate 

n > = -(ie) 2 E|G/A| 2 / 0^Tr{S>(P)S>(K-P)} (28) 

where the sum is over all the fermion species in the plasma. As usual we consider 
only SM fermions to be thermal. Using (jUJ) and ((7|) and some algebra this becomes 

IP = ~\G fA \ 2 -m 2 ^\j dp J- [(1 - n a )(0_ - n 2 )Q- + ^(0+ ~ "2)6+] . 

(29) 

To avoid clutter we have defined 

hx/2 = n(Ey 2 ) 
9 ± = 9(k ±E 1 ) 



Ei = \Jp 2 + mj 

E 2 = ^J(k-p) 2 + m). 

and we have used the identity 

e(po)(l-/(po)) = 0(po)-n(po) (30) 

The symbols 
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± 



1, if - 1 < (K 2 ± 2k Q E 1 )/2kp < 1 
0, otherwise 



impose kinematic constraints. 

It appears that IP diverges as k — > 0. But, in this limit, G + vanishes completely 
and the region of integration determined by 0_ collapses to a point given by K 2 = 
2k Ei. This gives a finite result. To evaluate the limit we put the A on-shell and set 
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K = (m, ). This is a good approximation since Y A jm A is typically on the order of 
0.05. Writing the infinitesimal region of integration explicitely, we have 



lim n> 



-4a 



\G 



,2 / m A 



I rq+2k 



71 



p 



ni){0. 



n 2 



(31) 



Where q = ^Jm\/A — m 2 . Finally, evaluating the integral, we obtain 



limhP 



-2am A 1 



2m 2 / 
m\ 



\ 



Am 2 
m\ 



(32) 



The imaginary part of the real-time self energy at zero external momentum is then 



ImlT 



1 (l _ e -/J*o) n > 



am A \G 



fA\ 



2m 2 / 
m 2 A 



i 



Am 2 
m\ 



1 



[m,A — 2m 



(33) 

There are several points worth noting. First, the self energy has no imaginary part 
for m A less than 2m j. This is expected for the following reason: at zero external 
momentum the processes that contribute to the boson's decay width are A — ► //. 
The additional processes A°f — > / and A°f — > f, corresponding to Landau damping, 
are kinematically disallowed. Thus, a minimum mass of m A > 2m f is required for a 
decay to take place. Instead of a power-law dependence on T, the thermal A width 
is supressed by a fermionic statistical factor evaluated at half the A mass. We could 
have deduced this fact as well by noting, again, that the only relevant process is the 
one involving the production of a fermion-antifermion pair and its inverse . These two 
processes should be accompanied, respectively, by statistical factors (1 — n(m A /2)) 2 
and n{mA/2) 2 with a relative minus sign between them. The total statistical factor 
would then be 

This is exactly the temperature dependence we obtained in (|33|). In Fig. Qwe plot 
the thermal corrections to the A width as a function of temperature. It is clear that 
at freeze-out (T ~ 5 GeV), the thermal effects are completely negligible. 



6 Conclusions 



In this paper we have identified several mechanisms which lead to heightened sen- 
sitivity of the WIMP annihilation cross section, and thus the relic density, to the 
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Figure 7: Thermal correction to the width of the A boson. 

values of the underlying theory parameters. We then investigated the consequences 
of finite temperature effects in regions of parameter space where these mechanisms 
are important. We first considered the issue of thermal corrections to parameters 
to which the relic density is exponentially sensitive, specifically the LSP-NLSP mass 
splitting in the stau coannihilation region. We found a thermal correction of about 
10" 4 to the mass splitting and a similar correction to the relative Boltzmann factor 
between the stau and the LSP at freeze-out. This correction is far to small to have 
any measurable consquenses for current or upcoming relic density observations. In 
the bulk region we considered the issue of corrections to the standard model fermion 
masses and the corresponding shifts in helicity suppressed annihilation cross sections. 
While the branching ratio for these processes is small, the fermions are light and 
can thus receive significant mass corrections. It is the interplay of these two effects 
which we investigated. We found that, for a typical point in the mSUGRA bulk, the 
thermal effects correct the fermionic mass operators by at most a factor of 1CT 4 at 
freeze-out. Finally, we considered thermal corrections to the width of the A boson 
in the A-funnel region. We found the value of this parameter to be exponentially 
suppressed by the ratio of the A mass to the freeze-out temperature and we provided 
a physical interpretation for this result. 

Although in the specific family of models which we studied all thermal effects 
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are considerably below the detection threshold of any current or planned observation 
of the relic density, we believe that generically these issues should be considered. 
For example, in any coannihilation scenario an exponential dependence on a small 
mass splitting is expected. Similarly, any resonant annihilation will exhibit hightened 
sensitivity to a small width associated with the virtual particle. Small corrections 
to these parameters cannot a priori be neglected. However, as we have seen, close 
analysis will often indicate that these effects are very small. 
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